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Abstract 



We show that the decay B° s — > J/ipf (9S0), / (980) — > 7r + 7r~ can be used to measure the CP 
violating phase in _B S mixing, —2/3 s , and estimate the sensitivity as ±0.050 rad, for 2 ftr 1 of 
LHCb data. After adding in the related B® — > J/ipt]', rf — > /ry mode, the sensitivity improves 
to ±0.044 rad. Use of these CP eigenstates obviates the need for a transversity analysis that 
must be used in the case of — > J r /i/xf) decays. 



1 Introduction 

While CP violation in B° decays has been measured unequivocally, attempts at determining the 
CP violating phase in B s mixing, — 2j3 s [1] by the CDF and DO experiments using B® — > J/ip(p 
decays [2], have given values much larger than Standard Model (SM) predictions, but with large 
enough errors that the measurements are not statistically significant, even when combined [3]. 
Since this mode is not a CP eigenstate, but one involving two vector particles an angular 
analysis is necessary [4]. These analyses, however, did not allow for the possibility of there 
being an S-wave component in the <ft mass region, thus possibly biasing the result and surely 
underestimating the error [5]. 

Since physics beyond the SM can or even should contribute virtual particles that interfere in 
the B s mixing loop, it is important to measure the mixing phase as precisely and in an unbiased 
manner as possible. It has long been known that the modes B® — > J/iprj^ > are CP eigenstates, 
and thus angular analysis is not needed. All of the 77 or rf decay modes contain at least one 
photon, whose reconstruction is much less efficient than charged particles in LHCb. One study 
has estimated the sensitivity using rf — > p7 [6], but the sensitivity is of the order of a factor of 
two worse than that expected using J/ipcj) mode. Other modes have also been considered [7]. 

In this note we describe the event selection, the backgrounds, and make an estimate of the 
measurement sensitivity of —2f3 si using a heretofore not considered mode B® — ► J/ipf (980), 
where /o(980) — ► 7r + 7r~. The dominant Feynman diagrams for the <fi and fo processes are 
shown in Fig. 1. There also are possible small contributions from penguin diagrams [8] and 
W-exchange [9]. 
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Figure 1: Decay diagrams for B s — ► J/ip<p, and B s — > J-0 Jo- 
There have been several estimates of the relative widths of these decay modes. Based on 
studies of hadronic decays Stone and Zhang [5] estimate that the relative width is 



There also is a non-/ 7r + 7r~ S-wave component that is very wide in mass, and contributes an 
additional 5% of the J j ifxft Ta ^ e using a narrow selection of ±90 MeV around the /q mass that 
was used in the sensitivity estimate [5]. 

Recently, the CLEO collaboration [10] has reported another estimate of Rf/$ using mea- 
surements of the semileptonic widths of — > foe + v and — > 0e + z/ at the endpoint of the 
four-momentum transfer range where the phase space is maximum, as suggested by Stone and 
Zhang. These measurements yield 



jr(D+-+/ (980)e+i/, f -> tt+tt- 



A note on notation, since it is clumsy to always refer to — 2(3 s we set this equal to 0/ to 
indicate that this is the CP violating phase measured in B s — > J/ipfo decays, which should to 
an excellent approximation be the same phase as measured in B s — ► J/ip<f> decays even in the 
presence of new physics in the mixing amplitude. 



2 Signal Selection and Optimization 

The present study is done using a Monte Carlo simulation of the signal and specific backgrounds. 
For the signal, a B s is produced in a 14 TeV proton proton collision conforming to the theoretical 
fragmentation function. It is allowed to decay and the quasi-stable particles then traverse the 
detector where they are subject to magnetic fields, multiple scattering in material, decays and 
hadronic interactions [11]. To optimize the selection, we generate a signal sample, and several 
background samples. We use a generic pp — ► J/ip X sample to predict the backgrounds when 
we have a J/ip present in the event, though these events should be eliminated since the J/ip is 
produced at the primary vertex (PV). To study the background from b decays we use a generic 
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b decay sample that also includes b — > J/ip X decays. However, some individual modes present 
specific problems, so we also generate them separately. Table 1 lists the different decay channels 
that have been generated, the number of analyzed events passing the geometrical cuts, and the 
geometrical cut efficiencies. 

We first generate the events and pre-select the ones within the geometrical acceptance: for 
exclusive b decays we insist that the charged B candidate decay tracks are below 400 mrad and 
larger than 10 mrad with respect to the beam line. In the case of the inclusive J/ip sample, we 
insist only that both muon tracks pass the above mentioned criteria. For the inclusive bb, we 
only require that the one of the two B mesons be pointed within 400 mrad of the beam line. 



Table 1: MC samples used in this study. 



Decay 


Number of events 


Geometrical cut efficiency 


B» - J/^(^)/o(7T7r) 

B° d - J/4>(w)K* (Kn) 
B+ -> J/^fi)K+ 

B° s -> J/i){w)i{f>l) 
B + ^ JMw)X 
B° d -> J/i>{w)X 
B° s - J/^{nn)X 
inclusive J /if) (////) 
inclusive bb 


60,565 
2,994,516 
1,928,159 
63,319 
2,831,431 
2, 924, 276 
720, 109 
4, 237, 046 
20,176,844 


16.4% 
17.3% 
17.9% 
16.4% 
20.4% 
20.4% 
20.2% 
19.7% 
43.7% 



Table 2 lists the production cross sections at 14 TeV predicted by Pythia/EvtGen, used in 
our Monte Carlo simulation. 1 

Table 2: Production cross sections at 14 TeV predicted by Pythia/EvtGen, used in our Monte 
Carlo simulation. The double arrow (i.e. =>■) indicates the existence of possible intermediate 
states, e.g. cr pp ^j/^ x includes (r pp ^ b ^j/i>x- The "prompt J/tp" cross section is cr prompt j/^ = 
v PP ~^j/ipx ~ & P p^b^j/ipx =0.266 mb. 



cross section 


mb 




102.9 


@ PP — >bb 


0.698 


<y pp ^j/%i}X 


0.286 


<y pp ^b^j/ipx 


0.0204 



To reconstruct the J/ipfo {J/fp — > , fo — > 7r + 7r _ ) candidates, we first insist that we have 

two opposite sign muon candidates that form a J/ip candidate and satisfy the "pre-selection" 
criteria listed in Table 3. Then we require two additional charged tracks, that satisfy other 
loose-pre-selection cuts, and are consistent with making a vertex with the two muons. The 
very loose pre-selection criteria are chosen with the aim to remove as little signal as possible, 
and reject large parts of the combinational background. This pre-selection is also applied to 
inclusive bb and inclusive J/ip events. 

1 We kept the LHC energy at 14 TeV in order to compare with other simulations involving <pf. For 10 TeV 
running, the final signal yields should be scaled by a factor of w 5/7. 
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Table 3: Cut values for the B° s — > J/ipfo pre-selection and selection. Also shown are the 
selection efficiency of each cut for the signal (absolute) and background (with respect to "pre- 
selection.") ACij is the likelihood ratio of species % relative to species j. The product of the pr 
values of J/ip daughters is used to help eliminate any residual minimum bias background. 



Cuts on the muons 


Pre-Selection 


Selection 


signal efficiency(%) 


BG cfliciency(%) 


A£ M7r 


> -5 


> -5 


97.2 




xLck/nDof 




< 5 


98.9 


85.9 




Cuts on the J/ip 


Product pt of daughters 




> 500 2 MeV 2 


100 


100 


Mass Window 


±42 MeV 


±42 MeV 


97.9 






Cuts on the pions 




> -10 


> -10 


99.2 






> -10 


> -10 


99.9 




xLck/ nDof 




< 4 


95.3 


65.2 


Min IPS 




> 3 


74.8 


1.5 




Cuts on the / (980) 


Sum pt of daughters 




> 900 MeV 


98.6 


57.6 


Mass Window 


±500 MeV 


±90 MeV 








Cuts on the B° s 


Max IPS 




< 5 


99.7 


67.5 


vertex fit x 2 


< 50 


< 17 


88.1 


51.8 


cos9 p 




> 0.99993 


80.4 


4.8 


Mass Window 


±300 MeV 


±50 MeV 


98.5 





2.1 J/fp — > selection 

Muon candidates are selected by requiring the track has hits in the Muon Chambers and 
satisfies the identification criterion that requires global A£ M7r > —5, and has Xtrack P er # °f de- 
grees of freedom (nDOF) of the Kalman fit of the track < 5 (Fig. 2 (d)). Specifically A£ P1P2 
is the difference in the log of the likelihood between the hypothesis that this particle is of type 
of p\ rather than p 2 (equivalent to the likelihood ratio). These cuts were studied in Ref. [12] 
using minimum bias events and are aimed at rejecting hadrons misidentified as muons due to 
random combinations of spurious hits in the Muon Chambers. The distributions of transverse 
momentum (p T ) for /i + vs. the minimum impact parameter significance (IPS) with re- 
spect to each primary vertex, 2 x 2 rack /nDOF and momentum of muon candidates are shown in 

2 The Impact Parameter Significance is defined as the track impact parameter with respect to a primary 
vertex divided by its error: IPS = ( ^2-). At the luminosity of 2 x 1CT 32 cm 2 s 1 about 35% of the events have 
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Fig. 2(a-e) from the signal and inclusive bb MC after "pre-selection." 

We then combine two opposite sign muon candidates to form a J/ip candidate. Fig. 3 
shows the signal and bb distributions for IPS, the flight distance significance (the distance from 
PV to reconstruction vertex divided by its error), the vertex fit x 2 an d invariant mass of J/ip 
candidates from the signal and inclusive bb MC after "pre-selection." We retain those events 
having |m MA1 — mj/^\ < 42 MeV corresponding to about a ±3cr interval (Fig. 3 (d)). These 
selections result in a relatively clean J/ip sample. 

2.2 / (980) — > 7T+7T- selection 

To select pion candidates we veto tracks identified as kaons with ACxn > 10 or muons with 
AC^ n > 10. The veto has > 99% efficiency for signal. Fig. 4 shows the distributions of for 
7r + vs. 7r~, the minimum IPS with respect to each primary vertex, x 2 rack /nDOF and momentum 
of pion candidates from the signal and inclusive bb MC after "pre-selection." 

We then combine two opposite sign pion candidates to form a / (980) candidate. Fig. 5 
shows the minimum IPS, flight significance, and vertex \ 2 distributions for /o(980) candidates 
from the signal and inclusive bb samples after "pre-selection." 

2.3 B® reconstruction 

A B° s is reconstructed by combining J/ip and /o(980) candidates. We require the four tracks 
from J /if) and /o(980) to be consistent with coming from one vertex as evaluated by examining 
the vertex fit x 2 , where the nDOF is 5. We also calculate the cosine of the angle between the 
B Q S candidates reconstructed momentum and direction from the primary vertex to B° s vertex 
(cos# p ), Fig. 6 shows the distributions of B° s candidates from the signal and inclusive bb MC 
after "pre-selection." 

2.4 Optimization Criteria 

To study how the <pf sensitivity depends on the signal statistics (S) and the background to 
signal ratio (B/S), we use a fast stimulation or "Toy" Monte Carlo. 

We generate and fit 400 "experiments" as a function of the B/S ratio, using <pf set to its 
predicted SM value of -0.0368. For each B/S point we fit the resulting <f>f distribution to a 
Gaussian. (This is the same fitter as used in section 6.) The r.m.s. width (a) of the Gaussian 
is taken as the 0j sensitivity. For our first point we fix B/S at zero and change only the 
amount of signal. Fig. 7 shows the error of 0/ (o^/) as a function of signal yield. The curve 
follows the expected behavior: a^ s oc 1/ \J~S. To understand the relation between o$ f and 
B/S, we fix S and only change B/S. Fig. 8 shows the the error of 0/ (c^) as a function 
of B/S, where the experiment contains signal and a long-lived background, with a lifetime 
fixed at 0.96 ps, as obtained from B — > J/ipX Monte Carlo. We fit the data with a shape 

a ( f >f oc \JT+~a x B/S. The value of a is 0.63±0.10, consistent with the ratio of background 
lifetime to B s lifetime, 0.64. For an exponential proper time distribution the lifetime can be 
quickly estimated using the mean of the distribution. For convenience we take a as the ratio 

more than one primary vertex, and we choose the smallest IPS. 
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(c) Minimum IPS of muons 
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Figure 2: The distributions of the signal (black solid) and bb background (red dashed) for muon 
candidates from J/ip. The curves in (a) and (b) show the product of Pt(^ + ) ■ Pt(h>~) = 500 2 
MeV 2 , used to eliminate minimum bias events in the preselection. The final cut in (d) is 
indicated by vertical blue dotted lines. 
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(a) J/ip minimum IPS (b) J/if> flight significance 
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(c) J /if) vertex \ 2 (d) J /if) mass 

Figure 3: (a)-(c): The distributions of signal (black solid) and bb background (red dashed) for 
J/ip candidates. (Numbers of events in the B and S distributions are equal.) (d) The J /if) 
mass distribution for signal, fit to a single Gaussian, with the selection cut indicated by vertical 
lines. 
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Figure 4: The distributions of signal (black solid) and bb background (red dashed) for pion 
candidates from f . The final cuts in (a)-(d) are indicated by blue dotted lines. (The number 
of events in the B and S distributions are equal.) 
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Figure 5: The distributions of signal (black solid) and bb background (red dashed) for /o(980) 
candidates. (The number of events in the B and S distributions are equal.) 
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(a) vertex fit \ 2 



(b) cos 6^ 




(c) IPS (d) Flight significance 



Figure 6: The distributions of the signal (black solid) and bb background (red dashed) for B® 
candidates. Final cuts on (a)-(c) are indicated by vertical blue dotted lines. 
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of mean values of the proper-time distributions of background and signal. (We have also found 
that the cf)f uncertainty due to the prompt J /if} background follows this formula.) 




v 

12 3 

S 



Figure 7: The error of 0/ as a function of signal number S. (S equal to one corresponds to 
about 2 fb _1 data.) The curve is the result of a fit to a^ f oc 1/yS. 



2.5 Selection Optimization 

Our goal is to maximize the 0/ sensitivity. Specifically, the selection cuts are chosen to max- 
imum ^g=g B where S (B) is the expected signal (background) number within a ±50 MeV 
mass window around the B° s nominal mass [13], and a is the ratio of the mean values of the 
proper-time distributions between the background and the signal. The cuts are based on a set 
of variables that show a marked difference between signal versus background. They include: 
(1) sum of absolute p T of f daughters (Fig. 4 (a,b)); (2) impact parameter significance of pions 
(Fig. 4 (c)); (3) track fit x t 2 rac k/ nD0F of pions (Fig. 4 (d)); (4) vertex fit x 2 of (Fig. 
6 (a)); and (5) the cosine of the angle between the B° s candidates reconstructed momentum and 
its flight direction (cos 6^) (Fig. 6 (b)). (We have not explored optimizing the selection based 
on a neural network or similar procedure as the tuning of such a method must be based on real 
data, and here we are trying to understand the sensitivity to first order. If other variables also 
have some discrimination power between signal and background, we will include them in the 
future.) 

In what follows, we normalize all the yields to 2 fb _1 data. We use the expected bb production 
cross-section a b i = 500 fib and o" promp t j/^ = 265.6 /ib as predicted by Pythia, and listed in 
Table 2. Since the Monte Carlo used 698 fib as the bb production cross-section, we scale the 
number of background events from bb to the value corresponding to a h i = 500 fib. To have more 
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Figure 8: The error of 0/ as a function of background to signal ratio B/S for a fixed S. The 

S corresponds to about 2 fb _1 data. The curve shows a fit to a^ f = AJl + a x B/S. Here 
a = P2 and A = PI. 

statistics in the background estimation, we enlarge the fo window from the nominal 90 MeV 
to 500 MeV, and the B® window from 50 MeV to 300 MeV. If the background is distributed 
linearly in the larger windows, we can scale the background number in the larger windows down 
by a factor of 33. The factor a may change among cut variables, so we obtain a from the signal 
and background proper-time distributions at each set of cut points. (We regard this as a minor 
point, and first analyses are likely to use a fixed value of a.) 

We realize that the lower mass region contains a large B° — > J/ipK*° sample when the 
kaon is misidentified as a pion (see Fig. 11a); there is also a small B® — > J/ip7r + 7T~ background 
(see Fig. 9), so we remove this background from the calculation of B as we will consider this 
background separately. The final cuts, called "Selection," are shown in the third column of 
Table 3. 

A selection cut of 900 MeV on Pt(tt + ) + pt(h~) serves to eliminate 42.4% of of the back- 
ground, costing only 1.4% of the signal efficiency 

3 Signal efficiency and event yields 

Table 4 shows the efficiencies computed for B° s — > J/ipfo, where: 

• e geo is the efficiency of the acceptance cut at generator level on the B° s daughters between 
10 to 400 mrad; 
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• e dct/ g co is the efficiency of that all B® decay products in the event are reconstructible (has 
enough MC hits for long track reconstruction) ; 

• ^rec/dct is the efficiency that the reconstructible events are actually reconstructed. 

• e se i/ rcc is the efficiency that the reconstructed events are actually selected. 

• e_L0) chlti, and €hlt2 are respective efficiencies for the LO, HLT1 and HLT2 trigger 
algorithms. 

Table 4: Efficiency for B° s -> J/^/„(980). LO, HLT1 and HLT2 refer to the three trigger steps. 



e geo [%] e dct / gco [%] 


^rec/det [%] 


Csel/rec [%] e L0 [%] 


CHLTl [%] 


CHLT2 [%] 


e to t [%] 


16.4 38.7 


71.4 


45.4 94.0 


95.0 


95.0 


1.75 



The overall trigger efficiency is 85%, for signal events with all decay products inside the 
detector. This is a high efficiency, especially compared with fully hadronic events, which are 
typically at least a factor of two lower. The total efficiency including HLT2 is the product of the 
individual efficiencies listed above, which is 1.75%. The yield for a decay channel is computed 

as 

Y = L- mt xa bi x2x f Bs x B vis x e tot , (3) 

where L int = 2 fb _1 is the assumed integrated luminosity, a h i = 500 /A> is the expected bb 
production cross-section, f Bs = (11.0 ± 1.2)% is the probability for a 6-quark to form a B s 
meson, and £> V1S is the total visible branching fraction, the product of all the individual branching 
fractions involved in the decay chain. (The factor of 2 arises because we can use either the b or b 
quark.) Assuming that B{B° S — > J/ipfo) x &(fo —> ii + 7i~) is quarter of B(B® — > J/ip<j)) x B(<f> — >■ 
K + K~), we have B V1S (B® — > J/ip (////) /o( 7r7r )) — 6.8 x 10~ 6 . After applying the efficiencies, we 
expect 26.1k signal events after HLT2. Table 5 shows the branching and 6-hadron production 
fractions assumed in the calculation of the yields and background levels discussed below. 



Table 5: Branching and Miadron production fractions assumed in the calculation of the yields 
and background levels. 



Branching fraction 


value 


estimated from 


Ref. 


B™(B» - J/^)f (nn)) 


6.8 x 10~ 6 




[5] 


£vi S(jB o ^ Jl^)rf{pri)) 


10 x 10" 6 


B» ^ J/i>K» 


[14] 




(5.25 ±0.24) x 10~ 5 




PDG08 [13] 


B™{B+ a - J/^{w)K+) 


(5.9 ±0.2) x 10~ 5 




PDG08 [13] 


fu 


(39.9 ± 1.2)% 




PDG08 [13] 


fd 


(39.9 ± 1.2)% 




PDG08 [13] 


fs 


(11.0 ± 1.2)% 




PDG08 [13] 
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4 Background Sources 



4.1 Prompt J/ift 

We use Monte Carlo corresponding to a 0.00125 fb _1 data sample and find 69 candidates passing 
the selection which are from prompt J/tp events in the enlarged fo and mass regions, before 
applying the trigger selections. From prompt J/ip events we expect B/S = (11 ± 1)%, where 
the error reflects only the uncertainty due to Monte Carlo statistics. 

4.2 Background from bb 

Our Monte Carlo simulation used B(b — > J/ipX) = lA6% (see Table 2), compared with the PDG 
value of (1.16 ±0.10)% [13]. Using the Monte Carlo value, we expect B/S = (26 ±2)% from bb 
background, where the error is only the statistical uncertainty from the Monte Carlo. This is 
reduced to (20 ±2)% using the PDG value. For the purpose of this note we will use the larger 
value, as that was what was used in previous studies of J/ip<p. We discuss specific background 
sources below. 

4.2.1 B u , d , s -> J/^X 

Fig. 9 shows the /i + /i~7r + 7r - invariant mass distribution from B u ^ s — > J/ipX. Fig. 10 shows 
the m nn distribution in the B® mass signal region. 

4.2.2 B° d J/tf>K*° 

When the kaon is misidentified as a pion, B® — > J/%pK*° decay can fake a B s ^ J/ip-K + Tc~ event. 
Fig. 11 shows the invariant mass distributions for B® and fo candidates from B® — > J/if)K*° 
decays. If the mass resolution does not deteriorate, we will see almost no events from this 
source; quantitatively, in 2 fb _1 of data, we expect < 10 events at 90% confidence level (C.L.) 
remaining in the fo and B® narrow signal regions. 

4.2.3 B+ J/ipK+ 

The .B+ — > J/%l)K + decay can fake a J/^7r + 7r _ when the K + is misidentified as pion and a 
random rr~ is combined. Fig. 11 shows the invariant mass distributions for B° s and f from 
5+ — > J/il)K + decays. We expect (1.29 ± 0.08) x 10 3 events contribute to the background 
before the trigger in 2 fb _1 data. 

4.2.4 B° s J/if)r]' 

We have studied rj' — > i]7i + 7i^ and rj' — > p7 and found only the latter could contribute as 
background. It is a dangerous background, because B® — > J/iprj' has opposite CP to the signal. 
Fig. 13 shows the invariant mass distributions for B® and fo from B° s — > J/ipi]', r( p7 
decays. Assuming B(B° S -> J/W) = | &{B° -> J/4>K°) [14], we expect (0.75 ± 0.07) x 10 3 
events contribute to the background, before the trigger, in 2 fb _1 data. Most of these, however, 
are not in the narrow signal mass ranges. 
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Figure 9: /i + yU~7r + 7r~ invariant mass distributions from B u ^,s J/ipX MC sample in (a) ±90 
MeV and (b) ±500 MeV windows around the /o(980) nominal mass. The expected signal 
distributions are superimposed. The non-resonant — > J/ip-K+n^ branching fraction in this 
MC generation is set to 2 x 10" 4 , about 30% of that of B° s -> J/ip(j>, <p -> K + K~. 
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Figure 10: The m n7T distribution in B° s mass signal region (\m^^ n7T — rriBo \ < 50 MeV) from the 
sample including all B decays with J/ip candidates, except signal. 
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Figure 11: Distributions from B® — > J/ipK*° decays when the kaon is misidentified as a pion. 
The vertical lines show the narrow mass windows. 
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Figure 12: Distributions from B+ — > J/ipK + decays where the K + is misidentified as pion and 
a random tt~ is combined. The vertical lines show the narrow mass windows. 

We also intend to reconstruct these events when we can find the photon from the rf — ► p7 
decay to use in a separate measurement of 0/ (see Sec. 6.4). We have not vetoed these events 
in our J/ipfo selection, in this study. 
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Figure 13: Distributions from B° s J/ipr]', 7/ — > p7 decays. The vertical lines show the narrow 
mass windows. 



4.3 Background Summary 

Table 6 shows the sources of the background and expected yields from 2 fb _1 data. The 
branching fractions and 6-hadron production fractions used for the calculation are shown in 
Table 5. The largest background arises from a J/if> combined with random tracks identified as 
pions, other than — > J/ipK*° and B® — > J/iprj'. These backgrounds (mass shape and proper 
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Table 6: Summary of background sources. The numbers quoted are calculated in the narrow 
mass windows — my ( 980 ) | < 90 MeV and \m^ n7T — m B o \ < 50 MeV. We assume the trigger 
efficiencies are the same for the signal and backgrounds. The bb background does not include 
the specific B° s1 £> + , B® channels listed here. 



Sources 


Yield from 2ft) 1 before trigger (xlO 3 ) 


B/S 


signal 


29.3 




background 






prompt J /if) 


3.3 ±0.4 


(11 ±1)% 


bb 


6.6 ±0.5 


(21 ± 2)% 


B° s -> J/W 


0.75 ±0.07 


(2.4 ±0.2)% 


B+ -> J/i)K+ 


1.29 ±0.08 


(4.2 ±0.3)% 


B° d -> J/ipK* 


< 10 @90%C.L. 




Total 


11.3 ±0.6 


(39 ± 3)% 



time distribution) can be estimated by reconstructing J/ip with two same sign-charged pions; 
we cannot use the B° s mass sidebands due to reflections from specific final states. We show in 
Fig. 14 the J/^tttt mass distribution for like sign and opposite sign pion pairs. In the B° s signal 
region, and even above, the distributions are in agreement. 




5100 5200 5300 5400 5500 5600 

M B> (MeV) 



Figure 14: Invariant mass distributions from from B° s — * J / ipmr candidates. The Monte Carlo 
sample consists of any B decay into a final state containing a J/if> where the modes B® — > J/ipf , 
B° s — > J/ipi]' and B® — > J/i/j7t + it~ have been explicitly removed. The histogram is for 7r + 7r~ 
pairs, while the crosses represent 7r ± 7r ± pairs. The vertical lines define the B® signal window . 

Furthermore, as a check, we processed a 5.6 million mini-bias event sample that satisfied the 
L0 trigger. Only 2 events passed the off-line selection, and they both contain real J/i[) decays. 
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5 Properties of selected signal events 
5.1 B® mass resolution 

Fig. 15 shows the reconstructed B® mass distribution for the signal. Performing a double- 
Gaussian fit gives a(M B o) = 19 MeV, where the average width a = \J(l — /2V1 + iWf , <Jx and 
cr 2 are the width of the core and second Gaussian, and f2 is the fraction of the second Gaussian 
[15]. 
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Figure 15: Reconstructed B° s mass for the signal. The HLT1 trigger has been applied. 



5.2 Proper time 

The B° s proper time is defined here as: 



where m is the reconstructed invariant mass, p the momentum and d the distance of flight 
vector of the candidate B® from the primary to the secondary vertices. If more than one 
primary vertex is found, the one that corresponds to the smallest IP of the B° s candidate is 
chosen. 

Fig. 16 shows the distribution of the difference of the reconstructed (t rec ) and MC true (t MC ) 
proper time. The average proper time resolution obtained from a double-Gaussian fit is 34 fs. 

Fig. 17 shows the proper time error estimate and the proper time pull trcc ~ f — . The estimate 
of the error on the reconstructed B® lifetime has a mean value of 27 fs and a most probable 
value of 25 fs. 

The time-dependent selection efficiency is shown in Fig. 18 for events before and after the 
trigger requirements. It is parameterized by the acceptance function defined as: 
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Figure 16: The distribution of the difference between the reconstructed (t rec ) and MC true 
(t ) proper time. The HLT1 trigger has been applied. 




Figure 17: Distributions of (a) the proper time error estimate and (b) the proper time pull. 
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where C is the selection efficiency at large decay time, and a and n are two parameters which 
govern the proper time dependence at small decay times. 




I (ps) t (ps) 



(a) untriggered (b) triggered 



Figure 18: The acceptance function as function of proper time. 



6 Sensitivity of — ► J/x/jfo Compared With J/^4> 

In order to perform this task a fast (Toy) Monte Carlo simulation program has been used that 
is based on RooFit. We simulate 400 LHCb "experiments" with 2 fb _1 per sample. For each 
experiment the program generates events taking into account the results obtained from the full 
MC simulation and performs a fit according to a likelihood function which includes the B Q S —B S 
mixing variables. The sensitivity to 0/ is taken as the root-mean-square width, a, resulting 
from a Gaussian fit to the <pf values from each simulation run. 



6.1 Likelihood function 

The time- dependent decay rates for initially pure B° s or B° s states in b — > ccs or b — > ccs 
quark-level transitions are given by the following expressions for decay into CP eigenstates: 

AT J . ATJ 



T(Bg — > fcp) oc e Fat jcosh — ^ r/jcos0/sinh — ^— + rjf sin (f)f sin(Am s t)| 

T(B° S — > fcp) oc e~ Tst jcosh — ^ ?7jcos0/sinh — ^ rjf sin <fif sin(Am s t)| (6) 



where t is the proper time, 77/ is the CP eigenvalue of the state fcp, and Ar s is the lifetime 
difference between CP+ and CP- eigenstates. Direct CP violation is neglected. Note that for 
J/tpfo: Vf = - 1 - 

The events are used to maximize a likelihood function (£) which is given by 



C= HP(m t ,tr,q t ), (7) 
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with 

p{m h tr, ft) = u g p^(m)pr(tr, ?o + a - /s ig )^ ks (m,)p t bks (tr c ), (8) 

where: 

• P,^ g (mi) and P^ g (rrii) are the probability density functions (PDFs) describing the de- 
pendence on reconstructed mass mi for signal and background events; 

• P t slg (t- ec , ft) is the PDF used to describe the signal decay rates for the decay time t] 00 , 
which depends on the tagging result at time t=0, ft (q = +1 if the signal meson is tagged 
as B®, q — — 1 if it is tagged as B° s , or q = if no tagging information); 

• P t bkg (tJ ec ) is the PDF describing the background decay rates, which do not depend on the 
tagging performance. 

• /sig is the fraction of the signal in the fitting region. 

The likelihood function includes distinctive contributions from the signal and the back- 
ground. For both, the PDF is a production of PDFs which model the invariant mass distribu- 
tion and the time-dependent decay rates. The PDF used for generating the B® mass spectrum 
consists of a double- Gaussian for the signal and a linear function for the background. The 
PDF for the proper time is modeled by a true decay function smeared by time resolution then 
multiplied by the decay time acceptance function. From Eq. 6, the true time decay function 
for the signal can be expressed as: 

R(t, q) oc e~ Tat jcosh — ^— + cos0/ sinh — ^~ qD sin0j sin(Am s t)| . (9) 

The effect of the wrong-tag probability w tag is included in the dilution factor D = (1— 2w tag ). We 
take this value as fixed from studies of other modes such as B° s — > Djn + , and use w tag = 0.334 

as our value. We take the wrong tag fractions as the same for B° s and B° s . 



6.2 Input to the fast Monte Carlo simulation 

The PDFs and their input parameters are shown in Table 7. The parameters of the signal are 
obtained from the distributions shown in Fig. 15, 16 and 18 (b). For the parameters of physics 
and tagging performance, we use the same as those used in B® — > J/i/xj) study [8], as shown in 
Table 8. The PDF of the signal proper-time is shown in Fig. 19 for both 0/ = —0.0368 (SM 
value) and 0/ = —0.2 rad, where the oscillations are apparent. 

The background's proper time distribution is obtained from B UjdjS — > J/ipX shown in Fig. 20. 
The average proper time is 1.93 ps for the signal and 1.20 ps for the bb background. We use 
B/S = 0.26 for background from b decays and B/S = 0.05 for that from prompt J/ip. 3 

3 The tagging efficiency for the prompt J/ip background is lower than the signal and bb background [8]. 
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proper time (ps) 



proper time (ps) 



(a) 4> f = -0.0368 rad (b) (p f = -0.2 rad 

Figure 19: The proper-time PDF of the signal for B° s (solid line) and T3° s (dashed line). 

6.3 Fit Results 

For this baseline study, we fit the Toy Monte Carlo only allowing cf)f and Ar to float, with the 
other parameters fixed to their input values. Fig. 21 shows the distributions of mass and the 
proper time for the events tagged as B° s and ~B° S from one experiment. The output of <pf from 
the fits is shown in Fig. 22 with signal only and background included. The sensitivity of 0/ 
is estimated to be (0.050 ± 0.002) rad. We also find that allowing the proper time acceptance 
parameters to float does not increase the uncertainty of </>/. 

6.4 Inclusion of B s — > J/i^r]' Events 

We have already discussed the the contribution of the J/tprj', rf — > p°7 events as a source of 
background in the J/ipfo sample. We also intend to reconstruct these events when we can find 
the photon from the rj' — > p°7 decay and add them into our final sample (with reversed CP to 
fo events). In our simulation we use all photons found in the electromagnetic calorimeter and 
those that convert in material in front of the magnet, provided that their pt is larger than 300 
MeV/c. The detection efficiency for photons in the solid angle of the detector is about 25%. 
In 2 fb _1 we estimate 5000 of such fully reconstructed events before the trigger. A previous 
analysis of this mode [6] concluded that an error in the measurement of 0/ of ±0.8 rad could 
be made with a 2 fb _1 sample. Adding the two modes together would give an error in 0y of ± 
0.044 rad. 

7 Systematic Errors 

We have studied several sources of systematic error. Recall the outputs of our fit nominally are 
4>f, AT and the time acceptance parameters (a and n). We actually determine, however, the 
product of the dilution D times 0y and use the value of D determined from other measurements. 
Thus the systematic error on D is fully correlated with the systematic error on <pf. 
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Table 7: The PDFs for the invariant mass and proper time describing the signal and background. 



Signal Double Gaussian (2G) 
2G(m;m ,5 m ,a 1 ,a 2 , / 2 ) 
m = 5368.4 MeV 
5 m = -22 MeV 
(Ji = 16.8 MeV 
a 2 = 33.2 MeV 
h = 0.09 



[R(t MC , q) ® 2G(f 



■ c -t MC ; / x,a*,a*,/ 2 t )]-e t (t rcc ;a,n) 
/x = 0.0008 ps 



a. 



0.0236 ps 
0.047 ps 
= 0.36 
a = 1.68 
n = 1.87 



TMrr 



Background 
from b 



First-order polynomial 



/i = 
(T* = 0.039 ps 



T 



bkg 



a 
n 



0.96 ps 
= 4.4 
= 2.9 



MC 



Background 
from prompt 

J/1> 



First-order polynomial 



2G(t rec - 1 



• e t {t vec ;a,n) 

H = 



/I 



0.11 ps 
2.4 ps 
= 0.17 
a = 4.4 
n = 2.9 



Now we will estimate the systematic error on uj tag . This parameter can be measured using 
a combination of other modes. One simple approach is to measure u tag using B° s — > D+n~ . We 
expect differences in the value of ui tag here and in the J / ipf mode because of different triggering 
in the hadronic and dimuon channels [16]. This difference can be estimated by simulation and 
checked using other modes. For example, we can use B° — > J/ipK*°, K*° — > K + ir~ to separately 
measure opposite side tagging. As input to a first estimate we decided to see the difference in 
Monte Carlo between the J/V'/o and D+ir~ B s final states. 

Our simulation of B Q S — > J/ipfo yields u tag = (0.326 ± 0.003), where the error is purely 
statistical. Using exactly the same tagging code on the mode B° s — > Dfn~, which of course we 
can and will use to measure u tag , gives a value of (0.309 ± 0.003), again the error is statistical. 
(We believe that we can measure Ut ag with excellent statistical precision, so we will ignore the 
statistical error.) The difference is 5.3% in ujt ag , or 10.6% in D. We argue that to first order 
we can use the Monte Carlo simulation to make the 5.3% correction to account for tagging 
differences between the two modes. Knowing that Monte Carlo is not perfect, we estimate that 
we will have « 25% error or ±1.3% uncertainty on uo tag after subjecting the Monte Carlo to 
various tests. Thus an systematic error of 2.6% on D should be achievable without superhuman 
efforts. We also expect that our belief in the Monte Carlo will change with time, hopefully 
improving. 

Two other sources of error have been investigated thus far using toy Monte Carlo. The 
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Parameters 


input values 


r s 


0.680ps _1 


A "P 


0.049 ps 


Aro s 


17.77PS- 1 


0/ = -2ft 


-0.0368 rad 


£ tag 


0.564 


^tag 


0.334 



Table 8: The input parameters of physics and tagging performance used for fast Monte Carlo 
simulation. 

general procedure is that we vary a parameter, or a shape parameterization, that has been 
fixed in the fit by plus and minus the expected error, and then repeat the fit. We record 
the difference of 0/ output between this fit and the nominal fit. Then the distribution of the 
difference is used to set the error by fitting to a single Gaussian. The mean of the Gaussian 
is taken as the systematic error for each particular source. To check if the systematic errors 
depends on the central value of 0/, we use both 0/ = —0.2 and 0/ = —0.736 for generation. We 
notice that the systematic error does depend on 0/ and cr^/0/ is approximately constant. We 
plan to measure the time resolution in real data from prompt J/xjj — > where we add two 

other tracks from the primary vertex consistent with the /o mass. We have seen that the pull 
distribution of time resolution for the signal and the prompt events, so defined, are identical. 
Our results are listed in Table 9. 

Table 9: Systematic error sources on 0/ 



Source 


Value 


<W0/ (%) 


Time resolution 


34 fs varied by ±5% 


2.5 


Misstag rate 


0.334±0.004 


2.6 


B s lifetime 


1.470±0.027 fs 


3.0 


Total 




4.7 



These errors are of comparable sizes. They will improve with increased statistics especially 
in the control channels. None of the systematic errors changes our estimation of statistics 
necessary for a significant effect, but merely scales the value of 0/. We are continuing these 
preliminary studies. 

Other more global sources of systematic errors due to B s production and B s decay rates are 
discussed in the Conclusions. 

8 Conclusions 

We expect 26,100 B s — > J/ipfo, /o — > 7r + 7r~ signal events in 2 fb _1 of accumulated LHCb data. 
Based on branching fraction predictions of resonant B® — > J / if) f (980) , /o(980) — > 7r + 7r~ and 
non-resonant B Q S — > J/ipn + -n~ , where the two pions are S-wave [5], we predict an error on the 
measurement of the CP violating parameter <pj of ±0.050 rad. Adding in the J/iprj', r( — > p7 
final state that we need to measure to estimate backgrounds, reduces the error on 0j to ±0.044 
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Figure 20: The background proper time distributions from (a) from B u ^ s — > J/ipX and (b) 
prompt J/ijj. 

rad. This is larger than the estimate using the J/ip<fi final state of ±0.03 rad, however the 
latter estimate does not consider the effect of a K + K~ S-wave. Initial indications are that 
taking the S-wave into account will increase the error by less than 15% [17]. Use of the J/ip 
plus scalar or pseudoscalar CP eigenstates removes the need for a complicated angular analysis 
and should provide, at minimum, a crucial check on the vector- vector result. Both methods 
have a large ~ ±25% systematic uncertainty on the predicted sensitivity due to uncertainties 
in the production cross-section and B s branching ratios. In addition, the estimate based here 
assumes that relative yield of the 7r + 7r~ S-wave in the fo mass region is 25% that of K + K~ in 
the cf) mass region. This introduces a considerable uncertainty of about ±40% in the relative 
rates and ±20% in the sensitivity, giving an overall uncertainty in of ±32% in the error in (pf. 
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Bs mass (MeV) proper time (ps) 



Figure 21: (a) Invariant mass and (b) proper-time distributions for the events tagged as B®, 
and (c) invariant mass and (d) proper-time distributions for those tagged as ~B^ S . The points 
with error bars are fast simulation data and the solid lines are fit functions. The dashed lines 
show the background contribution. 
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Figure 22: The 0/ output from fitting to 400 toy MC. 

9 Acknowledgments 

We thank the U. S. National Science Foundation for support. We also thank our LHCb col- 
leagues for useful discussions and review of this paper. 



28 



References 

[1] Note that S ^ -2ft, since -2ft = axg(V tb V*) 2 / (V cb V* s ) 2 , whereas S is arg(M 12 /r 12 ), 
with arg(M 12 ) = V tb V t * s ) 2 / (V cb V* s ) 2 , and arg(T 12 ) is a linear combination of (V cb V* s ) 2 , 
V cb V* s V ub V* s , and (KftV^) 2 . See A. Lenz, Nucl. Phys. Proc. Suppl. 177-178:81-86, 2008, 
arXiv:0705.3802v2 [hep-ph]. 

[2] D. Tonelli (CDF), arXiv:0810.3229 [hep-ex]; T. Aaltonen et al. (CDF), Phys. Rev. Lett. 
100, 161802 (2008), arXiv:0712.2397 [hep-ex]. V.M. Abazov et al. (DO), Phys. Rev. Lett. 
101241801 (2008) arXiv:0802.2255 [hep-ex] and J. Ellison (DO), arXiv:0810.1041 [hep-ex]. 

[3] G. Punzi, "Flavour physics at the Tevatron," presented at "The 2009 Europhysics Confer- 
ence on High Energy Physics," Krakow, Poland, July, 2009, to be included in the pro- 
ceedings. For older results see J. Charles (CKM Fitter) http : //www. slac . Stanford, 
edu/xorg/ckmf itter/TalkBuf f er/StatusOf TheCkmMatrixCapri2008.pdf, see also O. 
Deschamps (CKM Fitter) arXiv:0810.3139vl [hep-ph]. 

[4] I. Dunietz, H. R. Quinn, A. Snyder, W. Toki and H. J. Lipkin, Phys. Rev. D 43, 2193 
(1991). 

[5] S. Stone and L. Zhang, Phys. Rev. D79, 074024 (2009) [arXiv:0812.2832]. 

[6] D. Volyanskyy, "Study of B° s — > J/^rj'," presented at working group July, 2007, 
http : //indico . cern . ch/getFile . py/access?contribld=l&resld=0&materialld= 
slides&conf Id=9609. 

[7] Other such modes previously considered include B° s — > J/iprf( r y r y), B® — > J/ipr](ir + iT~iT ), 
B° s -> r] c (Ah)(f), B° s -> D+Dj L. Fernandez, Acta Phys. Pol B 38, 931 (2007) [LHCb-2006- 
047], B° s -> J/^V(P°7) D. Volyansky, and J. van Tilburg [CERN-LHCB-2007-027], and 
B° s -> J/H(W + ^) S- Jimnez-Otero [CERN-THESIS-2007-051]. 

[8] J. Albrecht et al., "Road map for the measurement of mixing induced CP violation in 
B° s -> at LHCb," LHCB-ROADM AP3-00 1 . 

[9] There is the possibility of an additional decay diagram involving VF-exchange, where we 
end up with a uu system with the J/ ip being formed from additional gluons. If this diagram 
is significant and if the fo has a relatively large uu piece in its wave function, it could be 
a source of an additional phase. Current estimates of the light quark content in the fo 
wavefunction, however, are small, see Ref. [5] and references therein. 

[10] K. M. Ecklund et al. (CLEO), "Study of the semileptonic decay D+ -> f (980)e + u and 
implications for B s -> J/^/o," [arXiv:0907.3201]. 

[11] For LHCb afncinados: we use DaVinci v20r3, based on DC06 Monte Carlo. 

[12] M. Calvi, et al. Lifetime unbiased selection of B° s — > J/ipcj) and related control channels: 
B° -> J/4>K*° and B + -> J/rj)K+ n , LHCb-2009-025, CERN-LHCb-2009-025 (2009). 

[13] C. Amsler, et al. (Particle Data Group), Physics Letters B 667, 1 (2008). 



29 



[14] K. Anikeev et al, "B Physics at the Tevatron: Run II and Beyond", arXiv:hep-ph/0201071. 

[15] We have not used a constraint on J /tp mass to improve the mass resolution because of the 
possibilites of introducing biases. This may be implemented in the future. 

[16] See M. Calvi, O. Leroy and M. Musy, "Flavour Tagging Algorithms and Performances in 
LHCb," LHCb-2007-058 ; CERN-LHCb-2007-058 and references therein. 

[17] P. Clarke, G. Cowan, F. Muheim, Y. Xie, " Determination of 2/3 s in B° s -> J/^jK+R- 
Decays in the Presence of a KK S-Wave Contribution," [arXiv:0908.3627]. 



30 



